Can Relic Superhorizon Inhomogeneities be Responsible for Large-Scale CMB Anomalies? 
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We investigate the effects of the presence of rehc classical superhorizon inhomogeneities during inflation. 
This superhorizon inhomogeneity appears as a gradient locally and picks out a preferred direction. Quan- 
tum fluctuations on this slightly inhomogeneous background are generally statistical anisotropic. We find a 
quadrupole modification to the ordinary isotropic spectrum. Moreover, this deviation from statistical isotropy 
is scale-dependent, with a ~ — factor. This result implies that the statistical anisotropy mainly appears on 
large scales, while the spectrum on small scales remains highly isotropic. Moreover, due to this — factor, 
the power on large scales is suppressed. Thus, our model can simultaneously explain the observed anisotropic 
alignments of the low-l multipoles and their low power. 



I. INTRODUCTION 

The current observation data support the standard 
ACDM model fl]]. Recently, however, there is an grow- 
ing interest in analyzing possible large-scale anomalies 
of CMB, from both theoretical and observational sides 

It appears that the lowest CMB multipoles are anoma- 
lous in two seemingly distinct aspects. Firstly, the angu- 
lar power Ci at the lowest £ is abnormally suppressed yfl. 
Secondly, they have an improbable directionality revealed 
by the fact that for a certain preferred orientation, one m- 
mode absorbs most of the power, which may imply the 
presence of an "axis of evil" flUt]. 

The statistical properties of perturbations carry the same 
symmetries as the background on which they are gen- 
erated. In the standard scenario, quantum fluctuations 
are assumed to be generated on a spatially homogeneous 
and isotropic background. Background spatial inhomo- 
geneities are neglected in most of the analysis. Indeed, 
a long period of inflation pulls all non-smooth classical 
initial conditions out of the horizon. However, if infla- 
tion lasts the minimal number 60 e-folds or so, perturba- 
tions with comoving wave-numbers of cosmological inter- 
est cross the horizon just during the earliest several e-folds 
of inflation. Thus, the relic inhomogeneities may leave 
some marks in the primordial quantum fluctuations. 

In this paper, we investigate the possible effects of a relic 
classical superhorizon inhomogeneity during inflation, es- 
pecially its effects on the statistical properties of quantum 
fluctuations. In a single superhorizon perturbation 

mode were considered to explain the power asymmetry of 
CMB on large scales. While in this work, we consider 
relic superhorizon inhomogeneities as background. Per- 
turbation theory in the presence of a spatially inhomoge- 
neous inflaton background value has been investigated by 
several authors before lfTlll . lfl2|] . Quantum fluctuations on 
this slightly inhomogeneous background are generally sta- 
tistical anisotropic. 

The deviations from statistical isotropy can take on 
many forms, which may correspond to different physical 



origins. One simple form was presented in 101, ^'(fc) = 

P{k) 1 + g{k){k ■ n)'^ . In this work, we find that 

the (leading-order) correction to the statistically isotropic 
power spectrum is of the ACW form but with a k- 
dependent factor (7(A;) ~ — Thus, the power spec- 
trum deviates from isotropy mainly on large scales (small 
fc), while remains highly isotropic on small scales. More- 
over, the spectrum itself is also suppressed on large scales. 
Thus, our model can simultaneously explain the observed 
anisotropic alignments of the low-£ multipoles and their 
low power 



n. MODEL AND BACKGROUND 

In this work we investigate general single scalar field 
inflationary models, with an action of the form S 



-g[R/2 + P{X, (/))], where X 



Now we consider a slightly inhomogeneous inflaton back- 
ground <j){t, x). Firstly, we assume the deviation from ho- 
mogeneity is very small, i.e., if (f>{t) is a spatial average of 
(f>{t, x) (e.g. in one Hubble volume), then we assume that 



^{t,x)-^{t) 



< 1. 



(2.1) 
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Secondly, we assume the inhomogeneities are superhori- 
zon. In other words, the typical comoving scale of these 
background inhomogeneities / is much larger than today's 
comoving Hubble scale, I ^ {aoHo)~^. Actually, it has 
been known long before that inflation can occur in the pres- 
ence of superhorizon initial inhomogeneities |9l|. In other 
words, the inflaton field initially should be smooth up to 
physical scales larger than ~ H^^. Thus, today's obser- 
vational universe can indeed inflate from an initial small 
patch with classical inhomogeneities with typical comov- 
ing scale I ^ (aoi?o)~^- 

Superhorizon inhomogeneities locally look like a gradi- 
ent, (t>{t, x) = (f){t, Xf)) + {x — Xf)) ■ S/(f>{t, Xq) + • • • . In 
this work, we neglect higher-order derivatives and treat 
as approximately constant (while (j) itself is indeed slow- 
rolling, for instance). 
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III. LINEAR PERTURBATIONS 

We split (j>{t, x) into background and fluctuation config- 
urations. 



4>{t, x) = <j>{t, x) + S(p{t, x) . 



(3.1) 



In the following we denote di(j){t,x) = Ai, which we have 
assumed to be constant. 

As a first-step investigation, we neglect the metric per- 
turbations. Thus the calculation is straight forward since 
all we have to do is to expand the action for the scalar field 
directly. According to ( 13.1b . simple Taylor expansion of 
P{X, (p) around the background cj) up to second order of c 



gives -iE^'^ 9^509^ (5(?i- 
where we have defined 



P 



(3.2) 



Note that due to the non-vanishing background gradient 
dicj), S*^ is not proportional 5ij any more. The presence 
of the additional term which is proportional diipdjcf) in S'^ 
breaks spatial rotational invariance and will be responsible 
for the generation of statistical anisotropies of the scalar 
perturbations Scf). 

After introducing a new variable x) = 

aV— the second-order action for the scalar 
field perturbations can be written as 



'(27r)3 2 



2Ve7 
aH 



{A ■ k)u'_f,Uk 



a 



A. 7(^-fc)' 



Now we take the scale factor as a{rj) = —l/Hrj. Re- 
markably, in this simplest case, equation i3.5i has an ana- 
lytic solution, 

uUv) - ^X^J^ e'^^'^+i^ei^^'V^i-ckrjYUia, z^+l, z) , 

(3.6) 

in which U{a,iy + l,z) is the confluent hypergeometric 
function, with 



1 



ic^k"^ - ^'y{A ■ k) 
A\A ■ k\^/T+7^ 



(3.7) 



z = —i\A- fc| + e7^?7^ , 
X^\A- fc|V7+^+ \/e7(^ • k) . 

Due to the presence of the factor A ■ k, the modes gen- 
erally depend on k rather than k ~ \k\. Here the coeffi- 
cient in ( 13.6b is chosen in order to satisfy the Wronskian 
normalization condition u'i^{i])'Uk{v) ~ '"fc ('7)'"fc('7) = 
In the limit Ai 0, it can be verified that ( 13.6b re- 

duces to the well-known functional form Ukitj) — > 

^e~^'^'^^') y/~7]Hl^\—ckri) , which describes nothing 
but the normalized mode function of a massive scalar field 
in pure de Sitter spacetime. Thus our mode solution ( 13.61 ) 
generalizes the standard homogeneous and isotropic back- 
ground to the case of the presence of superhorizon back- 
ground inhomogeneities. 



(3.3) 

where rj is comoving time, and a prime represents d/drj, 
Ti = a' /a, and M is an effective mass term, e is defined as 

e = (jp /H'^. Here we define two dimensionless parameters ''v'ia)' ^^''^ as z — > (when /3 > 2), we get 



B. Anisotropic Power Spectrum 

When all modes of cosmological interest exit the Hub- 
ble scale, that is, in 77 ^ limit, since t/(a,/?, z) — > 
£(£-1) 1-/3 



P 



X 



H^P 



7 



XX 



(3.4) 



In this paper, we assume 7 > 0. Note that the canonical 
case corresponds to = 1 and 7 = 0. In deriving ( 13.3b . 
we use the approximation that J^^" is a function of only rj, 
that is, we neglect the spatial dependence of E^". 



Uk{v)=A{k)e'^^^'ih'^-i^^^^{-chjr-'' , 

(3.8) 

where we have defined an anisotropic deformation factor, 

, (3.9) 



A{k) = ^^"'''^ / 4z|A.fc|^7 + e72 



r(«) 



„2^2 



A. Equation of Motion for the Perturbations and Solutions 

The classical equation of motion for the mode function 
according to the second-order action ( 13.3b is 

// , 2i^j{A-k) , 
+ ^ 

a a'^H^ J 

(3.5) 

As it can be seen, an anisotropic dispersion relation arises, 
which will be responsible for the generation of a statisti- 
cally anisotropic power spectrum. 



which is responsible for the anisotropic deformation of the 
power spectrum on large scales. Here a, v are given in 
( 13.71 ). In the standard senarios, A{k) is just A{k) = 1. In 
our case, A depends on fc, not only on its amplitude k, but 
also on it direction, more precisely, on A • fc (see FIG.[T]). 

We would like to investigate the leading order effect of 
Ai = di(i>. Takeing the limit — > in ( 13.91 ) and using 
Stirling's formula, we get 



A{k) = 1 - i- 



'y/e^{A ■ k) 



7(A • k) 

d^k^ 



v{v + \){'^v- 1) 
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£7 



(3.10) 
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FIG. 1: Deformation factor |^(fe)p as function of fc = |fc| 
and 9, where cos 9 = A-k. Parameters are chosen as c = 2, 
A = 7 = 1, e = 0.01, u = 1.52. We plot In k from -3 to 4, 
and 9 from to 7r/2. 



It can be seen directly from ( |3.1Qb that the anisotropic de- 
formation factor A{k) reduces to 1 when setting Ai = 0, 
as expected. Moreover, for those modes with wavenum- 
bers k perpendicular to A, i.e. A ■ k ~ 0, it will be 
shown that A{k) = 1, that is, these modes do not feel the 
presence of A, and thus get no corrections. Modes with 
A ■ k Q will get corrections from A, but in the limit 
A — > (or more precisely, typical scale of k much larger 
than A fc > \A\), it will also be shown that A{k) 1, 
which means that small scales get smaller anisotropic cor- 
rections from A. 



From (ITTOl i. 



\A{k)f 



1 



j{A-k) 



1) 



1) 



-n 



(3.11) 

It is interesting to note that the leading order correction 
term are quadrupole, there is no dipole correction. Thus, 
our result gives the anisotropic spectrum of ACW form i4||, 
P{k) = P[k) [l + g[k)[fi ■ kf] with 



9{k) 



-n 



(3.12) 
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FIG. 2: Leading-order quadrupole terms of |^(fc)| 
(for a given direction of fc). The black dashed curve 
is the exact result, the red curve corresponds to the 
quadrupole term given in l [3T0t . Here we have cho- 
sen that fc is parrallel to A without loss of generality. 
Parameters are chosen as ^4 = 10""^, c = 2,ns = 
0.96, e = 0.01. We plot k from 0.5 x 10~^ to 1. 



perhorizon scales reads, 

fc3 



AL(fc) 



\A(k)\^ 



H_y /ck 

2^) v ^ y cP.x 

2 I An^\\'2 



(3.13) 



Ai4k)\A{k)\\ 



where |yl(fc)|2 is given in ( 13.1 It . The general anisotropic 
power spectrum ( 13.131 ) is function of k and 9. In our model, 
its shape is similar to that of |A(A;)|2 (see FIG.[T]i- 



C. Suppression of Power on the Largest Scales 

The observational data appear to suggest a spatial modu- 
lation in the CMB spectrum. Especially, a vanishing CMB 
temperature auto-correlationi on the largest scales has been 
investigated by several authors f^. As mentioned before, 
this can be naturally explained in our model, due to the 
suppression behavior of |A(fc)|, g{k) ~ —1/k'^, (see FIG. 
|2] and FIG. Ell. 



In our model, g{k) < and it has an apparent k- 
dependence, g{k) ^ — l/fc^. This result has two impli- 
cations. First, the deviation from statistical isotropy of the 
power spectrum is not constant on all scales (while in |01 
and other works it is often assumed that g{k) is a con- 
stant g,). Statistical anisotropy of Ce spectrum mainly ap- 
pears at the lowest £'s, while the spectrum at higher ts 
remains highly isotropic (this can be seen clearly from 
FIGdl where for smaller fc, |.4(A;)| is highly anisotropic, 
while for larger fc, |.4(fc)| remains flat). Second, for a 
given direction k, since g{k) ^ — l/fc^, the spectrum it- 
self is suppressed at large scales (see FIG. |2]i. This phe- 
nomenon explains the loss of power at the lowest ^'s in the 
Ce spectrum. 

Finally, the power spectrum of scalar perturbation on su- 



IV. CONCLUSION AND DISCUSSION 

In this work, we studied cosmological perturbations in 
the presence of a superhorizon inhomogeneous inflaton 
background. We find that the leading-order correction 
to the standard isotropic power spectrum is a quadrupole 
anisotropy. Moreover, the deviation from isotropy is gen- 
erally scale-dependent. More specifically, the quadrupole 
correction term has a suppression factor g{k) ~ — l/fc^. 
And thus our model can simultaneously explain the ob- 
served anisotropic alignment of the low-£ multipoles and 
their low power. If the indication of these large scale CMB 
anomalies were to be confirmed by future experiments, the 
work presented here would improve our understanding of 
inflation and the very early universe. 
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FIG. 3: Suppression of the power on the largest scales 
(for a given direction k). The dashed line is tradi- 
tional isotropic spectrum with Us — 0.96. From top 
tobottem, A = 0.3 x 10"^, 0.5 x 10"^, 10"^, 0.3 x 
10^^,10^^. Parameters are chosen as e = 0.01, 
J = H — 1, c = 1.2. The range of k is from 10"'^ to 
1. 

In our work, the inhomogeneous background affects the 
perturbation through the non-linear structure of P{X, (f) of 
X. A similar effect would arise when one consider back- 
reaction of superhorizon perturbations through the nonlin- 
earity of the Einstein equations ifioll . 

In our approximations, we neglected the backreaction 
of the inhomogeneous inflaton background to the geome- 
try, which has to be considered in a more rigorous analysis. 
In that case, a non-FRW metric may arise. For instance, in 



ifTsll . anisotropic inflatonary models and perturbations are 
investigated, where an anisotropic spectrum is also pro- 
duced. In a FRW background, the quantization procedure 
is well undersood for canonical variables on sub-Hubble 
scales, where they evolve adiabatically, and thus the adi- 
abatic vacuum is well-defined. In this work, we choose 
the mode solution (13.6b by simply noticing that ( I3.6l l re- 
duces to the standard mode solution in the Ai — > limit. 
It should be noted, however, that the frequency term in the 
second-line of (I3.5l l is divergent when 77 ^ — oo, due to the 
presence of the anisotropic term ^{A ■ fc)^/ {aH)^. Thus, 
the presence of a WKB regime and also the quantization 
procedure in this case are non-trivial in general (see e.g. 
discussions in Iil5il ). We would like to investigate these 
problems in future work. 

Indeed, the work presented here is a general framework, 
while the results of our model are generic. A more de- 
tailed study including metric perturbations and possible 
non-FRW background geometry is needed. We would like 
to present it in a companion work lfl3ll . 
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